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Abstract  The  compound-Gaussian  model  is  often  used  in  radar  signal 
processing  to  describe  the  heavy-tailed  clutter  distribution.  The  important 
problems  in  compound-Gaussian  clutter  modeling  are  choosing  the 
texture  distribution  and  estimating  its  parameters.  Many  texture  distribution 
models  have  been  proposed  [1],  [2]  and  their  parameters  were  typically 
estimated  using  the  (statistically  suboptimal)  method  of  moments,  see 
[2].  In  this  paper,  we  develop  maximum  likelihood  (ML)  methods  for  jointly 
estimating  target  and  clutter  parameters  in  compound-Gaussian  clutter.  In 
particular,  we  estimate  (i)  the  complex  target  amplitudes,  (ii)  covariance  matrix 
of  the  speckle  component,  and  (iii)  the  texture-distribution  parameters.  Several 
existing  texture  models  are  considered:  (i)  gamma  (leading  to  the  well-known  K 
clutter  distribution  [1],  [2]),  (ii)  lognormal,  and  (iii)  Weibull.  Motivated  by  the  robust 
regression  model  in  [3],  we  also  develop  a  complex  multivariate  t  distribution 
model  for  the  clutter.  We  utilize  the  expectation-maximization  (EM)  algorithm 
to  estimate  the  unknown  parameters.  Numerical  integration  is  typically  needed 
to  compute  the  conditional  expectations  in  the  expectation  (E)  step  of  the  EM 
algorithm;  here,  we  employ  the  Gauss  quadratures  to  perform  this  integration. 
Interestingly,  the  proposed  complex  multivariate  t  distribution  model  does  not 
require  numerical  integration,  allowing  for  remarkably  simple  estimation  and 
detection  algorithms.  We  will  also  compute  Cramer-Rao  bounds  (CRBs)  for 
the  unknown  parameters  and  demonstrate  the  performances  of  the  proposed 
methods  via  numerical  simulations. 
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ABSTRACT 

Compound-Gaussian  models  are  used  in  radar  signal  pro¬ 
cessing  to  describe  heavy-tailed  clutter  distributions.  Im¬ 
portant  problems  in  compound-Gaussian  clutter  modeling 
are:  choosing  the  texture  distribution  and  estimating  its 
parameters.  Many  texture  distributions  have  been  studied 
and  their  parameters  were  typically  estimated  using  sta¬ 
tistically  suboptimal  approaches.  We  develop  maximum 
likelihood  (ML)  methods  for  jointly  estimating  target  and 
clutter  parameters  in  compound-Gaussian  clutter  using 
radar  array  measurements.  In  particular,  we  estimate  (i) 
complex  target  amplitudes,  (ii)  spatial  covariance  matrix 
of  the  speckle  component,  and  (iii)  texture  distribution  pa¬ 
rameters.  We  consider  two  existing  texture  models,  log¬ 
normal  and  gamma,  and  propose  an  inverse-gamma  tex¬ 
ture  model,  leading  to  a  complex  multivariate  t  clutter  dis¬ 
tribution.  Parameter-expanded  expectation-maximization 
(PX-EM)  algorithms  are  developed  to  compute  the  ML  es¬ 
timates  of  the  unknown  parameters.  For  lognormal  and 
gamma  textures,  Gauss  quadratures  are  utilized  to  im¬ 
plement  the  estimation  algorithms,  whereas  the  inverse- 
gamma  texture  model  does  not  require  numerical  integra¬ 
tion,  thus  yielding  remarkably  simple  estimators.  We  study 
the  performance  of  the  proposed  methods  via  numerical 
simulations. 

1.  INTRODUCTION 

Compound-Gaussian  models  have  been  used  to  characterize 
heavy-tailed  clutter  distributions  in  radar  as  well  as  to  model 
speech  waveforms,  fast  fading  channels,  and  various  radio 
propagation  channel  disturbances,  see  [1]  and  references 
therein.  Important  problems  in  compound-Gaussian  clut¬ 
ter  modeling  are:  choosing  the  texture  distribution  and  esti¬ 
mating  its  parameters.  Many  texture  distributions  have  been 
studied  and  their  parameters  were  typically  estimated  using 

The  work  of  A.  Nehorai  and  J.  Wang  was  supported  by  the  Air  Force 
Office  of  Scientific  Research  Grant  F49620-02- 1-0339,  the  National  Sci¬ 
ence  Foundation  Grants  CCR-0 105334  and  CCR-0330342,  and  the  Office 
of  Naval  Research  Grant  N000 14-01- 1-0681. 


(statistically  suboptimal)  method  of  moments,  see  [2],  In 
this  paper,  we  present  maximum  likelihood  (ML)  methods 
for  estimating  target  and  clutter  parameters  in  compound- 
Gaussian  clutter. 

In  Section  2,  we  introduce  the  measurement  scenarios 
with  lognormal  [2],  gamma  [2]— [4],  and  inverse-gamma  tex¬ 
ture  models 1 .  For  these  three  models,  we  develop  parameter- 
expanded  expectation-maximization  (PX-EM)  algorithms  to 
compute  the  ML  estimates  of  the  unknown  parameters  (see 
Sections  3.1,  3.2,  and  3.3,  respectively)  and  evaluate  their 
performance  in  Section  4. 

2.  MEASUREMENT  MODEL 

We  extend  the  radar  array  measurement  model  in  [6]  to  ac¬ 
count  for  compound-Gaussian  clutter.  Assume  that  an  n- 
element  radar  array  receives  P  pulse  returns,  where  each 
pulse  provides  N  range-gate  samples.  We  collect  the  spatio- 
temporal  data  from  the  tth  range  gate  into  a  vector  y(t)  of 
size  to  =  nP  and  model  y(t)  as2  (see  [6]  and  [7]) 

y(t)  =  AX<f>(t)  +  e(t),  f=l,...,JV,  (1) 

where  A  is  an  m  x  r  spatio-temporal  steering  matrix  of 
the  targets,  #  =  [0(1),  0(2), . . . ,  0(iV)]  is  the  temporal 
response  matrix,  X  is  an  r  x  d  matrix  of  unknown  com¬ 
plex  amplitudes  of  the  targets,  and  e(t)  is  additive  noise. 
Here,  we  assume  that  the  additive  noise  vectors  e(t)  are  in¬ 
dependent,  identically  distributed  (i.i.d.)  and  come  from  a 
compound-Gaussian  probability  distribution,  see  e.g.  [1]- 
[4]  and  [8] — [10].  We  now  represent  the  above  measurement 
scenario  using  the  following  hierarchical  model:  y(t)  are 
conditionally  independent  random  vectors  with  probability 
density  functions  (pdfs): 

Py\u{yif)\u(t)-,X,S)  =  exp  |  —  [y(t)  —  AX<fi(t)]H 
•[w(f)X]_1-[j/(f)  -  AA'0(f)]|  j \nu(t)X!\,  (2) 

1  The  proposed  inverse-gamma  texture  model  leads  to  a  complex  multi¬ 
variate  t  clutter  distribution,  which  generalizes  the  multivariate  t  distribu¬ 
tion  in  e.g.  [5]  to  the  complex-data  measurement  scenario. 

2  A  special  case  of  the  model  (1)  for  rank-one  targets  (i.e.  scalar  X )  in 
compound-Gaussian  clutter  was  considered  in  [8]. 


where  “ H ”  denotes  the  Hermitian  (conjugate)  transpose,  £ 
is  the  (unknown)  covariance  matrix  of  the  speckle  compo¬ 
nent,  and  u(t),  t  =  1, 2, . . . ,  N  are  the  unobserved  texture 
components  (powers),  modeled  as  i.i.d.  non-negative  ran¬ 
dom  variables.  We  consider  the  following  texture  distribu¬ 
tions: 

•  lognormal:  In  u(t)  follow  a  Gaussian  distribution, 

•  gamma:  u(t)  follow  a  gamma  distribution,  and 

•  inverse  gamma:  l/u(t)  follow  a  gamma  distribution. 
Our  goal  is  to  compute  the  ML  estimates  of  the  complex 
amplitude  matrix  X ,  speckle  covariance  matrix  £ ,  and  tex¬ 
ture  distribution  parameters  from  the  measurements  y  = 
[y(  1)T,  y( 2)T , . . . ,  y(N)T]T .  In  the  following,  we  present 
parameter-expanded  expectation-maximization  (PX-EM)  al¬ 
gorithms  for  ML  estimation  of  these  parameters  under  the 
above  three  texture  models.  The  PX-EM  algorithms  share 
the  same  monotonic  convergence  properties  as  the  “clas¬ 
sical”  expectation-maximization  (EM)  algorithms,  see  [11, 
Theorem  1],  They  outperform  the  EM  algorithms  in  the 
global  rate  of  convergence,  see  [11,  Theorem  2], 

3.  ML  ESTIMATION 

3.1.  PX-EM  Algorithm  for  Lognormal  Texture 

Assume  that  the  unobserved  texture  component  follows  a 
lognormal  distribution  (see  also  [2]);  equivalently,  /3(f)  = 
lnu(f),  t  =  1,2, ...  ,N  are  Gaussian.  We  further  assume 
that  /3(f)  have  zero  mean  and  unknown  variance  a2\  hence, 
the  unknown  parameters  are  0  =  {X,  £,a^}.  [Note  that 
the  mean  of  /3(f)  can  be  chosen  arbitrarily  because  the  speckle 
covariance  matrix  £  is  unknown.] 

We  now  develop  a  PX-EM  algorithm  to  estimate  0  by 
treating  /3(f),  t  =  1, 2, . . . ,  N  as  the  unobserved  (or  miss¬ 
ing)  data  and  adding  an  auxiliary  (dummy)  parameter  pp 
[the  mean  of  /3(f)]  to  the  set  of  parameters  0.  Under  the 
expanded  parameter  model,  the  pdf  of  /3(f)  is 

Pp(P{t)-,  M/3,  <?})  =  ^  exp{-[/3(f)  - pp]2/{2a2p)}.  (3) 

We  also  define  the  augmented  (expanded)  set  of  parameters 
0a  =  {X,  £a,  ap,  pp},  where  £  =  exp(pp)  ■  £a.  The  PX- 


EM  algorithm  for  this  model  consists  of  iterating  between 
the  following  PX-E  and  PX-M  steps: 


PX-E  Step: 

Compute 

Ti(v,0®)  = 

1  N 
-  ■  5] 

t=l 

'E/3|j/[exp[— /3(f)]  y(f); 0a5]}, 

(4a) 

TMrf)  = 

i  N 

t=l 

‘E /3|v [ exp[  /3(f)]  y(f); 0a] ]|, 

(4b) 

T3(y,e^)  = 

i  N 

N  ’  51  {</>(*)</>(*)" 

*=  1 

'E/3|j/[exp[— /3(f)]  y(f); 0a] ]}, 

(4c) 

ii 

<£> 

i  N 

^■^Eply[p2(t)\y(ty,eW}, 

t=l 

(4d) 

h{y;0^)  = 

1  N 

-•^E/3|y[/3(f)|y(f);0W], 

t=  1 

(4e) 

where 

o(:] 

=  {X^,£^,{a})^,pf} 

(4f) 

is  the  estimate  of  0a  in  the  ith  iteration  and  (4a)-(4e)  are 

computed  using  (5)  (below). 

PX-M  Step: 

Compute 

x(i+1)  = 

[Ah(S«)-1A]-1Ah(5W)-1 

•T1(y,0«)T3(y,0«)-1, 

(6a) 

£(i+l)  = 

5(<)  +  [/„l-Q(i)(5(i))-1]T1(y,0i; 

0) 

■T3{y-,e\ 

T,  (6b) 

U{i+1)  - 

M/3 

t6{y;eP), 

(6c) 

(4)"*"  = 

f4(y,0ii))-(M?+1))2, 

(6d) 

2/(i+1)  = 

exp (m?+1))  •  ^+1)> 

(6e) 

E/3|y[exp[-/3(f)]|y(f);0a] 

^  Ez=i  hi  '  ^(M/3,cr§)_(m+1)  ■exp{-vi(pp,a20)-1  ■  [; y(f )  -  AXcf>(t)]H  £~1[y(t)  -  AX<t>{t)}} 

Ya=i  hi  ■  Viivt 3,cr§)_m  •exp{-t>z(^/3,cr^)-1  •  [y(f)  -  AXcj>(t)]H ££1[y(t)  -  AX<j>(t)]} 
E/3|y[/3fe(f)|y(f);  0a] 

^  Ez= 1  hi  •  (VZcrpxi  +  p0)k  ■  Vi(pp,a2p)-m  ■  exp{-vt(pp,  af}^1  ■  [y(f)  -  AXcj>{t)]H £~1[y(t)  -  AX<p(t)]} 
Ef=i  hi  ■  vi{pp,a2p)~m  ■  exp{—vi(yp,  a/g)-1  •  [y(f)  -  AX<p(t)\H X^lyit)  -  AX<f>(t)]} 

where  vi(pp,  cr|)  =  exp(s/2apxi  +  pp)  and  k  €  {1,2}. 


where 


PX-E  Step: 


Compute 


S{i)  =  T2(y,0«) 

-Ti(y,eW)  -niy-G^)-1  ■T1(y,0(j))H,  (6f) 

QW  =  A[AH  {S^)~l  A]-1  AH .  (6g) 

The  above  iteration  is  performed  until  XW,  rW,  and  (<r|)0) 
converge.  Here,  Im  denotes  the  identity  matrix  of  size  m, 
We  now  discuss  computing  the  conditional  expectations 
in  (5).  First,  recall  the  Gauss-Hermite  quadrature  formula 
[12,  Ch.  5.3]: 


L 

f(x)  ■  exp(— a;2)  dx^^hi 

i—i 


(7) 


where  fix)  is  an  arbitrary  real  function,  L  is  the  quadra¬ 
ture  order  (determining  approximation  accuracy),  and  xi 
and  hi,  l  =  1, 2, . . . ,  L  are  the  abscissas  and  weights  of  the 
Gauss-Hermite  quadrature  (respectively),  tabulated  in  e.g. 
[13],  Using  (7),  the  Bayes  rule,  equations  (2)  and  (3),  and 
change-of-variable  transformation  x  =  (/3  —  yp)/{\/2crp), 
we  obtain  the  approximate  expression  in  (5). 


3.2.  PX-EM  Algorithm  for  Gamma  Texture 

We  now  model  the  texture  components  u(t),  t  =  1,  2, . . . ,  N 
as  gamma  random  variables  with  mean  one  (as  in  e.g.  [4]) 
and  unknown  shape  parameter  v  >  0;  hence,  the  unknown 
parameters  are  6  =  { X ,  27,  //  }.  (The  shape  parameter  v  is 
also  known  as  the  Nakagami-m  parameter  in  the  commu¬ 
nications  literature,  see  e.g.  [14,  Ch.  2.2. 1.4].)  This  choice 
of  the  texture  distribution  leads  to  the  well-known  K  clutter 
model,  see  [2]  and  [4]  and  references  therein. 

We  develop  a  PX-EM  algorithm  to  estimate  0  by  treat¬ 
ing  u(t),  t  =  1,  2, . . . ,  TV  as  the  unobserved  data  and  adding 
an  auxiliary  parameter  fiu  [the  mean  of  w(f)]  to  the  set  of 
parameters  Q.  Under  this  expanded  model,  the  pdf  of  u(t) 
is  [for  u(t)  >  0] 


1  /  JS  \  ^ 

Pu(u(t);v,n)  =  =-^-f — )  u(f)ly_1exp  [-  vu{t)/yu\ 

(8) 

where  T(-)  denotes  the  gamma  function.  Hence,  the  aug¬ 
mented  parameter  set  is  0a  =  {X,  27a,  v,  yu},  where  27a 
and  27  are  related  as  follows:  27  =  yu  ■  27a.  The  PX-EM  al¬ 
gorithm  for  the  above  expanded  model  consists  of  iterating 
between  the  following  PX-E  and  PX-M  steps: 


Tiiy-,0^)  =  ^ 

TMrf)  =  ^ 
TMrf)  =  ^ 

*4  (rf)  =  ^ 
h(v,oP)  =  ^ 


N 

t=  1 

N 

t= l 

N 

t= 1 


I  y(t);eW],  (9a) 

I  (9b) 

(9c) 


•E«lw[uW  1  I 

N 

■  ^E„|j/[lnu(f)  |  y{t)\d%\  (9d) 

t=  1 
N 

•£Eu|w[u(t)|y(i);0«]t  (9e) 

t= 1 

where  0^  =  {XW,  27a’\  v^l\  is  the  estimate  of  0a 
in  the  zth  iteration  and  (9a)-(9e)  are  computed  using  (10) 
(below)  with  g(u(t))  =  u(f)_1,  In  u(t),  and  u(t). 

PX-M  Step:  Compute 


X(?;+1)  =  [Aff(S'(i))-1A]-1Aif(S'(i))-1 

•T1(y,0«)T3(y,0«)-1,  (11a) 

S(i+D  =  +  [Im-Q^iS^)-1}T1{y,d «) 

•T3(y;  efr'T^y,  0^)H  [Im  -  Q«(S«)-i]*  (llb) 
ti+1)  =  h(y,e (<)),  (11c) 

27^;+1^  =  ^;+1)  •  27ii+1),  (lid) 


where 

S{i)  =  T2(y,0 «) 

-T 1(y,0«).T3(y;0«)-1-T1(y,^)H,  (lie) 

q(0  =  a[Ah  (S(i>)_1  A]-1  Ah  ,  (Ilf) 

and  find  that  maximizes 

v(i+ 1 )  =  arg m.ax  |  —  lnT(tz)  +  v\n  v  —  vh\\t§{y,  0^)] 
+vU{y,6{i))  -  z/|.  (11c) 

The  above  iteration  is  performed  until  XW,  27 W,  and  i/W 
converge.  The  conditional-expectation  expression  (10)  is 
obtained  by  using  the  Bayes  rule,  equations  (2)  and  (8),  and 
change-of-variable  transformation  x  =  i m/y.  The  integrals 


/0°°  g{xyu/v)  -py|u(y(f)  |  xy,u/v\  X,  27a)  •  a:1"  1  exp(-a:)  dec 
/0X5p!/|«(2/(^)  I  xyu/v\  X,  27a)  •  a:1"”1  exp(-x)  da; 


E«ii/[ff(uW)  I  yW;0 a 


(10) 


in  the  numerator  and  denominator  of  (10)  are  efficiently  and 
accurately  evaluated  using  the  generalized  Gauss-Laguerre 
quadrature  formula  (see  [12,  Ch.  5.3]): 


f(x)-x1'  1  exp(— x)  dx 


L 


^2wi(is-l)  f(xi(v-l)), 

i=i 

(12) 


where  f{x)  is  an  arbitrary  real  function,  L  is  the  quadrature 
order,  and  xi(y  —  1)  and  wi(y  —  1),  l  =  1, 2, . . . ,  L  are  the 
abscissas  and  weights  of  the  generalized  Gauss-Laguerre 
quadrature  with  parameter  v  —  1. 

The  computation  of  i/t+i  1  requires  maximizing  (11c), 
which  is  performed  using  the  Newton-Raphson  method  (em¬ 
bedded  within  the  “outer”  EM  iteration,  similar  to  [15]). 


3.3.  PX-EM  Algorithm  for  Inverse  Gamma  Texture 

We  now  propose  a  complex  multivariate  t  distribution  model 
for  the  clutter  and  apply  it  to  the  measurement  scenario  in 
Section  2.  A  similar  clutter  model  was  briefly  discussed  in 
[10,  Sec.  IV.B.3],  where  it  was  also  referred  to  as  the  gener¬ 
alized  Cauchy  distribution.  Assume  that  w(t)  =  u(t)_1,  t  = 
1,  2, . . . ,  N  are  gamma  random  variables  with  mean  one  and 
unknown  shape  parameter  v  >  0.  Consequently,  u(t)  fol¬ 
lows  an  inverse  gamma  distribution  and  the  conditional  dis¬ 
tribution  of  y(t)  given  w(t)  is  py\u(y(t,)  |w(t)_1;  X,  X),  see 
also  (2).  Integrating  the  unobserved  data  w(t)  out,  we  ob¬ 
tain  a  closed-form  expression  for  the  marginal  pdf  of  y(t): 


py(y(t);X,£,v) 


T(v  +  m) 

\ttS\  •  r(t/)  •  vm 


■{ 


i  + 


}—v—m 

,  G3) 


which  is  the  complex  multivariate  t  distribution  with  loca¬ 
tion  vector  AX<fi(t),  scale  matrix  X,  and  shape  parameter 
v.  Here,  the  unknown  parameters  are  9  =  {X,  X,  v } .  We 
first  estimate  X  and  X  assuming  that  the  shape  parameter  v 
is  known  and  then  discuss  the  estimation  of  v. 

Known  u:  For  a  fixed  v,  we  derive  a  PX-EM  algorithm  for 
estimating  X  and  X  by  treating  w(t),  t  =  1, 2, . . . ,  ./V  as 
the  unobserved  data  and  adding  an  auxiliary  mean  parame¬ 
ter  for  w(t),  similar  to  the  lognormal  and  gamma  cases  dis¬ 
cussed  in  Sections  3.1  and  3.2.  Here,  the  resulting  PX-EM 
algorithm  consists  of  iterating  between  the  following  PX-E 
and  PX-M  steps: 

PX-E  Step:  Compute 

w^{t)  =  (v  +  m)  ■  |i/  +  [y(t)  —  AX^ 4>(f)]H 

•[  (14a) 


for  t  =  1, 2, . . . ,  N  and 


Til)  =  ^  (14b) 

1  N 

T 2W  =  (14c) 

t= 1 

1  N 

T3(i)  =  (14d) 


PX-M  Step:  Compute 

X(i+1)  =  [Aff(5(,;))-1A]-1AH(S'(i))"1T1(i)(T3(i))C1(15a) 
X(i+A  =  {s,w  +  [/TO-gW(s,W)-1]-T1(i)(T3w)-1 


i  N 

<T®)H  [Ira  -  Q{i)(S^)-Y}/[^ 


(15b) 


where 


Qd)  =  A[AH  (S^)-1  A]-1  Ah  . 


(15c) 

(15d) 


The  above  iteration  is  performed  until  X^  and  27  W  con¬ 
verge.  Denote  by  X(°°\i/)  and  X(°°\v)  the  estimates  of 
X  and  X  obtained  upon  convergence,  where  we  emphasize 
their  dependence  on  v. 

Unknown  v :  We  compute  the  ML  estimate  of  v  by  max¬ 
imizing  the  observed-data  log-likelihood  function  concen¬ 
trated  with  respect  to  X(i/)  and  X(v)\ 


N 


v  =  arg  max  ^  In py  (y (t) ;  X(oo)  (i/) ,  I7(oo)  (v) ,  v) , 


(16) 


see  also  (13). 


4.  SIMULATION  RESULTS 

The  numerical  example  presented  here  assess  the  estima¬ 
tion  accuracy  of  the  ML  estimates  of  X,  X,  and  the  shape 
parameters  of  the  texture  components.  We  consider  a  mea¬ 
surement  scenario  with  a  3-element  radar  array  and  P  =  3 
pulses,  implying  that  m  =  9.  We  selected  a  rank-one  target 
scenario  with  =  1,  t  =  1,2 ,N,  complex  target 
amplitude  X  =  0.207  •  exp(j7r/7),  and 

A  =  b{vu)  (g>  a(t?), 

where  b(w)  =  [1,  exp(j2nzu),  exp(j4-7r-n7)]T  with  normal¬ 
ized  Doppler  frequency  vj  =  0.42,  and  a(i9)  =  [1,  exp(j27r$), 
exp(j47n?)]T  with  spatial  frequency  d  =  0.926.  Here,  ® 
and  “T”  denote  the  Kronecker  product  and  transpose,  re¬ 
spectively.  The  speckle  covariance  matrix  X  was  generated 
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Figure  1:  MSEs  for  the  ML  estimates  of  X,  and  E  as 
functions  of  N  under  the  lognormal  texture  model. 

using  a  model  similar  to  that  in  [16,  Sec.  2.6]  with  1000 
patches;  the  diagonal  elements  of  E  were  10.17.  Our  per¬ 
formance  metric  is  the  mean-square  error  (MSE)  of  an  es¬ 
timator,  calculated  using  2500  independent  trials.  The  or¬ 
der  of  the  Gauss-Hermite  and  generalized  Gauss-Laguerre 
quadratures  was  L  =  20. 

We  first  study  the  performance  of  the  ML  estimation  of 
the  clutter  and  target  parameters  for  lognormal  texture  in 
Section  3.1.  Here,  the  variance  of  / 3{t )  was  set  to  cr^  =  0.5. 
Fig.  1  shows  the  MSEs  for  the  ML  estimates  of  X  and  <r| 
and  the  average  MSE  for  the  ML  estimates  of  the  speckle 
covariance  parameters3,  as  functions  of  the  number  of  range 
samples  (snapshots)  N. 

We  now  examine  the  performance  of  the  ML  estimation 

bn  particular,  we  average  the  MSEs  of  the  elements  of  Re{  vechf  E ) } 
and  T iri  (  vechf  L  )  }7  /  .  Here,  the  vech  and  vech  operators  create  a  single 
column  vector  by  stacking  elements  below  the  main  diagonal  columnwise; 
vech  includes  the  main  diagonal,  whereas  vech  omits  it. 
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Figure  2:  MSEs  for  the  ML  estimates  of  X,  v,  and  E  as 
functions  of  N  under  the  gamma  texture  model. 

for  gamma  texture  in  Section  3.2.  We  have  set  the  shape 
parameter  to  v  =  2.  Fig.  2  shows  the  MSEs  for  the  ML 
estimates  of  X  and  v  and  the  average  MSE  for  the  ML  esti¬ 
mates  of  the  speckle  covariance  parameters,  as  functions  of 
N. 

Finally,  we  show  the  performance  of  the  ML  estimation 
for  inverse  gamma  texture  in  Section  3.2.  Here,  the  shape 
parameter  was  set  to  v  =  4.  Fig.  3  shows  the  MSEs  for  the 
ML  estimates  of  X  and  v  and  the  average  MSE  for  the  ML 
estimates  of  the  speckle  covariance  parameters,  as  functions 
of  N. 
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Figure  3:  MSEs  for  the  ML  estimates  of  X,  v,  and  L1  under 
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Summary) 

Develop  maximum  likelihood  (ML)  methods  for  jointly 
estimating  target  and  clutter  parameters  in  compound- 
Gaussian  clutter  using  radar  array  measurements. 

We  estimate: 

(i)  complex  target  amplitudes, 

(ii)  spatial  covariance  matrix  of  the  speckle  component, 

(iii)  texture  distribution  parameters. 


